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二 维 非 线性 对 流 扩散 方程 的 特征 混合 
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摘 E 本 文 构造 了 求解 非 线 性 对 流 扩散 方程 的 两 重 网 客 算法 ， 该 算法 首先 是 在 步 长 为 HH 的 粗 网 格 上 求 


解 一 个 非 线 性 问题 ， 再 利用 粗 网 格 解 得 到 一 个 线性 问题 并 在 细 网 格 上 求解 一 个 线性 问题 。 理 论 分 
析 与 数值 计算 表明 ， 该 算法 不 仅 消除 了 数值 振荡 现象 ， 还 极 大 地 提高 了 非 线性 对 流 扩 散 方程 的 计 


算 效 率 。 
关键 词 : 对 流 扩散 方程 ， 特 征 混 合 有 限 元 ; 两 重 网 格 法 ， 收 敛 性 
分 类 号 : AMS(2000) 65N30 中 图 分 类 号 : О241.82 文献 标识 码 : A 
1 引言 


本 文 考虑 用 特征 混合 有 限 元 求解 非 线性 对 流 扩 散 方 程 
Du 


co +b-Vu-V-(a(u)Vu) = f(u), (хт the Q~ J, 
u(x,0) = u(x), TER, (1) 
u(x,t) = 0, (x,t) € OQ x J, 


其 中 Q CRARAC WHNA RK MON, MAT = (0,T], Т > 0. с = c(2,t), b= 
(ду (x,t), ba(z, t)), a(u) = a(u,z,t), Хи) = Хи, z,t)e 

对 于 p > 1, LQ) 为 定义 在 2 上 的 标准 Banach 空间 ， 且 具有 范 数 | pe WNO) HA 
Fi YE || - |. AY Sobolev 22 [8], 9 (-,-) HA LO) FM AB, W = LeO, V= H(div;Q) = 
{у є (12(90))?2: у-у є L2(Q)}. FFE (1) HEL FRE: 

(a) 0< c <celz,t) <c* < +оо, 0 <а, < а(и) < а < +оо; 


(b) [b/c] +] 7 -(b/c)| < Cis 
Of(s д? f(s 

© Se + ge о, 
Ж с, ct. an a С, 及 C2 都 是 正常 数 。 

方程 (1) 具有 广泛 的 工程 应 用 背景 h-9j， 它 可 以 描述 质量 、 热 量 的 输 运 过 程 及 反映 扩散 过 程 
等 众多 物理 现象 ， 构 造 稳定 、 快 速 实用 的 数值 方法 ， 对 于 实际 问题 的 解决 有 着 重要 的 理论 和 实 
用 价值 。 众 所 周知 ， 用 标准 有 限 元 和 有 限 差分 法 求解 对 流 占 优 扩散 方程 经 常会 产生 非 物理 震 
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荡 或 者 数值 扩散 " 引 ， 因 此 克服 这 一 现象 的 许多 数值 方法 便 应 运 而 生 。 混 合 有 限 元 法 最 初 是 为 
解决 固体 力学 问题 而 提出 B10， 现 在 已 应 用 于 流体 力学 和 弹性 力学 的 诸多 方面 。 该 方法 受到 普 
遍 重 视 ， 是 因为 不 仅 能 够 求 出 未 知 变量 的 数值 解 ， 同 时 还 能 求 出 未 知 变量 的 导数 值 n1343 。 两 
重 网 格 算法 是 近 几 年 才 新 出 现 的 求解 非 线性 方程 的 高 效 算法 。 文 [13] 前 述 了 如 何 通过 收敛 的 粗 
网 格 解 及 泰勒 展 式 获得 细 网 格 的 Galerkin 有 限 元 解 ， 文 [14-16] 已 将 此 方法 用 于 相关 问题 的 求 
解 ， 并 取得 理想 的 加 速 收敛 结果 。 

本 文 针 对 非 线 性 对 流 扩散 方程 构造 了 特征 混合 有 限 元 两 重 网 格 算法 ， 该 算法 除了 体现 特征 
线 法 和 混合 有 限 元 法 的 特点 以 外 ， 最 重要 的 是 可 以 加 快 非 线 性 对 流 扩散 方程 的 收敛 速度 ， 是 一 
种 稳定 、 高 效 的 数值 算法 。 


2 ”两 重 网 格 算法 的 构造 


A (x,t) = Vel, t) 十 1b(z,t)|2， 并 设 与 特征 方向 7 = r(x, t) 相伴 的 算 子 为 e(z, t) ёч +b. 
Vu, 其 中 
ð _ о с(шт)д | b(z,t) | 


r yeda was (2) 
于 是 ， 方 程 (1) 可 以 改写 为 特征 形式 


(x,t) —V-(a(u)Vu) = flu), (2,t)€ 2x J, 
u(x,t) = 0, (x,t) € 29 х J, (3) 


u(x, 0) uo(z), rE Q, 
并 给 出 (3) НЕВЕ 5502008 (u,T,®) EW xV xV 满足 


(v3, w) 十 (У 'Ф,ш) та (f(u), w), wew, 


(Г,у) =(u,V-v), ve, (4) 
(®, v) = (а(и)Г, у), veV, 
其 中 
T=-Vu, Ф=а(и)Г. (5) 


取 时 间 步 长 At 并 在 时 刻 妇 = nAt, n=0,1,---,N, At = Т/М 求解 (4)。 在 时 刻 t = t R 
征 导 数 通过 以 下 方式 离散 


u(x,t”) ~ ulz, t?) 


ди" Š 
(927) ЧИ) (ж — ®)2 + (At)? 
u(x,t”) — u(z, t?—1) =o u” — ar} 


с At Deo? 


(6) 
其 中 
bz, 刀 -) 


э Ck oa 
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这 里 可 能 会 出 现 xzEQ 的 情形 ， Te iF: R 一 60 是 z 到 69 的 投影 算 
Ў, даз ЖР F(z) 的 对 称 点 。 


же. | УЧ), 260, 
2р1(Е(2)) – fg), ZER. 


对 于 ze QQ， 车 EQ， 当 At 充分 小 时 ， 总 有 5 € QQ， 所 以 该 延 拓 总 是 有 意义 的 。 

对 区 域 呈 进行 直径 为 hh 的 拟 一 致 多 边 形 齐 分 得 到 有 限 维 空间 Wh с W 和 Vi с V， 记 这 样 
的 前 分 结果 为 A; 并 称 为 细 网 格 。 这 样 就 可 以 将 (4) 的 时 间 离 散 特 征 混合 有 限 元 解 定义 为 ， 对 
Fn=0,1,---,N, Ж (ut, TR, Фл) є (Wh, Va, Ул) 使 得 


(о wn) + (V+ OR, wa) = (FUR) ш), wn € Wi, 
(TR, уһ) = (uk, V va), vn Є Vh, (7) 


(Фр.уһ) = (alu) TR ур), Vn Є Vh, 


HP ug = urls, t”), Gf? = url, t!) 

为 了 加 快 非 线 性 问题 (7) 的 收敛 速度 ， 这 里 采用 两 重 网 格 算法 求解 。 首 先 采 用 非 线 性 迭代 
ЖН Ag C A, 上 的 特征 混合 有 限 元 解 骆 ，I 和 及 惠风， 然后 再 求解 一 个 线性 系统 的 
Ж (ap TR, Фл) € (Wh, Vn, Vn) 满足 


(cro wn) + (V E wn) = (Flug) + РО) (ар — ар) wa), шк Wi 


(Гу, va) = (GR, У-у), Vn € Vh, (8) 
(OP, уһ) = (aluh) Eg, ут) + (a (uh) Th (aR — uh), ур), Va € Vh, 
其 中 心 << Н. 


两 重 网 格 算法 只 是 将 类 似 于 牛顿 的 迭代 法 限定 在 粗 网 格 上 ， 然 后 在 需要 求解 的 细 网 格 上 只 
求解 线性 化 后 的 问题 (8) 即 可 ， 可 以 起 到 很 好 加 速 收敛 作用 。 


3 ”两 重 网 格 算法 的 误差 分 析 


在 进行 两 重 网 格 算法 的 收敛 性 分 析 时 ， 首 先 引 进 离散 空间 和 投影 算 子 ， 然 后 讨论 (7) 的 解 的 
收敛 性 ， 再 利用 (т) 的 结论 分 析 (8) 的 收敛 性 。 
对 (7) 的 特征 混合 有 限 元 解 进行 误差 分 析 ， 需 要 用 到 Raviart-Thomas 空间 的 如 下 假设 


М.уһЄ Wa, va E Va. (9) 
BTN nk, SFR LO HRMS" = 9(,t")， 引 入 范 数 
N 1/2 
lole = (X atier) . (10) 
n=} 


利用 L 范 数 定义 空间 Wi ЖП Va OH (07, Г) є Wi x УШК 
(u", wh) = (Uh, wh), wh E Wh, (11) 


(Г, vn) = (ÒR, Va), va € Vh. (12) 
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同时 还 需要 定义 投影 In : (H1(Q))? 一 Vi, WE 
(У Ф", ш") = (У-(П„Ф”),шһ), wh E€ Wh. (13) 
而 且 ， 投 影 还 满足 唯一 性 
(и",У-уһ) = (âR, Y Vn), Von € Vh. (14) 


再 设 有 限 元 空间 Wr AV, RAW TREIE: 对 于 ve HQ) Rv e HQ) x HQ) 


[у = Mavi) < Су, 1<r<k+1, (15) 
lu — Gal] < С.А", | – alloc < С.л", 1<r<k+1, (16) 
IE = Ên < CHEIR, [E — Falo < СІ, 1<r<k+1, (17) 


其 中 与 Wi 及 Vi 的 多 项 式 的 次 数 有 关 。 这 些 假设 对 于 Raviart-Thomas 空间 总 是 成 立 的 。 为 
分 析 方 便 我 们 还 假设 问题 (1) 的 解 是 充分 光滑 的 而 且 可 以 达到 上 十 1 阶 光滑 性 。 
HFE V, 和 Whp， 当 vp Є Ул, un E€ Wnr AB И 


[уло < valh, 1 < 11А". (18) 


щр= 2, oo0 及 1 <r<k 十 1(k > 1) 时 ,根据 标 准 有 限 元 的 误差 估计 还 有 [2 


[и — 2 ||(о,т:12) + Allu — @л||»(го,т;нъу < СА Jull »(о,т;н”)- (19) 
(и — ûn) (и — ûn) 
| ot L?(0,T;L?) T h| ot L?(0,T;H") 
т Ou 
= Ch (leeren |27 bun) 9 
3.1 混合 特征 有 限 元 的 误差 分 析 
5183.108) сє L®(0,T; L?) WA 
ar"? < + KA)’, 
其 中 下 为 非 负 常数 ， T 
b(x, t” 
4 = 42) = a(z 一 aay): 
在 时 刻 t= 1", FA (4) A (7) 得 到 
чу-е) 
HV: —\-Ф,шһ) = (f(u") — (иһ), wn), wa E Wa; (21) 


(Г? — T}, va) = (u” — up, V Va), Va E€ Vh, 


(Ф” — BR, vp) = (а(и”)Г” — а(ий)Гр,ур), Vh Є Vh- 
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Sne = u?’ — ar, н" = ur Е ar, Cn = 了 一 下 及 和 = oF – Пьф”, 通过 (12)-(14) 得 到 


(wn) + (V-x",wn) = (yur ~ rte wp) 
+( st w) + (F(u) = f (un), wn), wh Є Wh, (22) 
(C7, Vn) = (д, У : уһ), уһ Є Vhs 
(х®, уһ) = (Br — nB”, уһ) — (а(и")Г" — а(ир)Гй уһ), Va E Vh. 
在 任 一 点 ze Q， 由 中 值 定 理 可 知 存在 一 点 om" 使 得 flur) - /(иў) = f’'(a”)(u® — ир), ТЇЙ 
Ри") – Fu) = f(a ул" — /'(о”)и^*. (23) 
同 理 存 在 8" (т) 使 得 
аш)" — а(иў)Гў = а'(8")°Г" — а'(8")и”Г” + а(ик)(Г" — PR) — а(иһ)с". (24) 
(22) P, Sun = и", Vn =X", уһ = C” HAG (23), (24)， 得 到 
nb” шй Ё" st n nn pn 
(cn 2 —, nu) + (a(uR)e”, c") 
ди" 一 eee ы 
= (ит At м") + (e At a”) 
+P — f(a, ш") + (ПАФ — Ф", С") + (a (B°) r”, с") 
7 
(а (Gp, 6") + (afp E” — ÎR), c") = DT (25) 
i=l 
于 (25) 的 左 端 ， 由 引 理 3.1 有 
n пт-1 
(co —, p) + (аби), с") 
2 т-112\ү _ ck n-1))2 |2 
> © (рит? ри) — Sah IP + al? (26) 
利用 [14] 中 的 结论 及 引 理 3.1 AAA A tI Тү, Te 
„ди? pu” art 5 
[711 = I(v Or e At А *) <o| дт? EE [а І. (27) 


n 1—1 
[Ta = П a) | < (и +") + Сие + о. (28) 


由 条 件 (а)-(с) 及 Holder 和 Young 不 等 式 可 以 给 出 Ts — Ту 的 估计 如 下 


ITs} < C(u”? + 1902), (29) 
Ts| < СФ" — ПФ" |2 + eller |, (30) 
тъ < CULM т"? + eller? < Clin? + eller i?, (31) 
Te] < CUE Ыш"? + eller? < С? + eller, (32) 
т < CHE” — ÊRI? + єс". (33) 
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将 (26)-(33) LA (25), Фе < а, /8, С =min(c,,ax)， 用 2At 乘 以 两 端 并 将 nn 由 1 加 到 m(m < 
N)， 再 注意 到 Ho = 0， 推 得 


Па" + У (Acie |?) 


< (|| а? + > + lorlzAt+ > Inn Pat 
n=1 
+ > С° – 2де У 1” – пФ" |22). (34) 
n=1 


引用 离散 的 Gronwall 引 理 ， 再 由 (15)-(17) 及 上 式 得 到 
н” + 的 мс") < O(At + Atta + вен), (35) 


Sh=O(A)M, MAS AAS ји" ит < jv? 上 pz 及 ?一 TR S e+e -eR 
结合 (15)、(16)、(18) 得 到 


lu” ари (У дат" = rpe)” 


n=1 


m 1/2 та К 1/2 
< ll (У acter?) + ln + (У даг" - êr) 
n=1 n=1 
< C(At+ А®). (36) 
这 样 便 得 到 如 下 定理 。 


定理 3.1 如 果 条 件 (aj-(c) IZ, Hh =O/(At), WHF (7) EX (ит, Te, BP) є Wh x 
Va x Vas 必 存 在 与 h, At 无 关 的 С, 使 得 


ju” — ur] + [Г — Гл < C(At + Һ®), (37) 


HHi<m<N,k>1. 
3.2 混合 特征 有 限 元 两 重 网 格 算法 的 误差 分 析 
下 面 分 析 两 重 网 格 算法 的 收敛 性 ， 这 个 方法 分 两 步 进 行 : 
第 一 步 ， 在 粗 网 格 Ан “Ey 求解 如 下 非 线性 问题 解 得 (иҗ, Г, OF) € (Wu, Vi, Ун) 


(o BE wn) + (У Фр wa) = (fu) он), wh Є Wa, 
(Th va) = (ий,У-ун), Ун € Vi (38) 
(OF ‚ун) = (aluh) h VH), va Є VH. 


第 二 步 ， 在 细 网 格 Au 上 ， 求 解 线性 系统 (8) 得 到 两 重 网 格 解 (r „Гр, фр) є (Wh, Vas Va) 
定理 3.2 ”如 果 条 件 (a)-(c) з, Ал = O(At)， 则 对 于 由 (8) 定义 的 (аз, TR, BP) є Wi x 
Va x 记 ， 必 存在 与 hh 和信 t 无 关 的 C， 使 得 


ju” – aM || + IE -El < C(At+ А + H*-2), (39) 
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其 中 k>1,M<N。 
证 明 令 
6% = йр 0р, у= Гр Ё. O= Фа – ПФ", 


ъ=, Ш (4) ME (8), HO we =E, vn = 0", уһ = ү" AR 
САС 09 


+ (ort en) + (Fu) — Flu) РОДОВ — uh); 


(40) 
("б") = ("0-6"), 
(0", л") = (Ф" — ПьФ", ү") — (а(и")Г" — aluh) ER — a! (uh T (GR — ин), 7"). 
由 假设 (a)-(c) Ж Taylor 展 式 可 知 存在 a?， 使 得 
Fu") = Fuh) — Д — ий) 
= идут" РОДЕ + SF" (aR Ca" — uh). (41) 
同 理 ， 由 Taylor 展示 得 到 
au) = аби) + (uy (u" — uf) + УЭ (ym и) (42) 
于 是 
абат)" — aah PR — а (uh EE OR — ч) 
= атг" — ÈP) + a(uh (ER — PR) а)" — "УГ — T” + T”) 
жаи) ER — гд) + FO — ag?) Be. (43) 


3 (41). (43) 代入 (40) 得 到 


(EHET e) + (июл) 


n n _ pnl п _ 
= (êE е SE en) (оон) + (родот е), е) 
+5 (F"(GE Cu" а £") + (ПФ Ф", у”) + (а(и”)(Г" — Рр), 7") 
H(A (URI (ГИ Г") 7”) + (a! (us yal", э") = (a (URE — ЁЗ, ") 


= (a (uh) Ch — Г), у") — (а (ир) "Г", у") + (a (we) (u” -uh CR- Г), 7") 


13 
HG ~ wh) PE") = Ут: (44) 
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左 端 项 估计 同 (26) 得 到 


ee 


Cx We т 
> элу (IEP - Пе") - SE e~ + ate, 


一 到 的 估计 同 (27)-(28) 


mi < ofS At + ат? 


L2(tn-1 tr; L?) 


а < э E +) + SE pae + Hr. 


其 余 项 的 估计 如 下 
031 < C(ai? + IEP). 
Ж (37) PKA = Н, т = т, ЖЕН (15)-(18) 得 到 
[Za] < Си" — uf)? ||? + En? < CH-2(H* + At) + er |, 
|T5| < CI” — 11," |? + elly"|!?, 


ITs] < СГ” – TRIP + е, 


Т < Сїт" coll” — PEI? + ella? < CREDIT” — гу |? + lly IP, 


Tg] < CIES I? + Elly" ll? < Clin"? + ela”. 
BEN EJE], n= 1,2,.… ,NN 中 取 值 最 大 的 一 个 ， 于 是 由 (19) 及 Holder 不 等 式 得 到 
IT] < CI — ао < Che tla Wes — Ре 1877. 
再 由 (18) 得 


Tiol < Се" — PIE +"? < CIEI? + ет", 
Tia] < Сг" г"? + єй"? < Chery? + ella IP, 
[Tol < Clu" — аз) ER – ЁР)? + е 

< CH (H* + At)? (Га — Г"? + AAH) + ет, 


ITis| < Cr — а )2 + ell” l? < Сн-(н* + Ab4 十 clanll2. 


(45) 


(46) 


(47) 


(48) 


(49) 
(50) 
(51) 
(52) 


(53) 


(54) 


(55) 


(56) 


(57) 


(58) 


将 (45)-(58) 代入 (44)， 取 适当 的 s， 得 用 2At 乘 以 两 端 并 将 n 由 1 加 到 М < N， 再 注意 
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到 上 =0， 有 


м 
eM? + У (412) 


n=l 


м M M 
At)? + "[? + "At+ У 12де 
Ba 2-19 | у | D | 


5 of | SF 


M M 
+ Уг" – PRAtt+ УФ" – 1,8"? At + HH" 十 At 


n=1 n=l 


sD дц" г +H- цен дар" ih = Pall 


二 
+H~?(H" + At) af > дг — гч]? + +) К (59) 
对 (59) 中 右 端 最 后 三 项 再 做 估计 ， 得 到 
Wem ||? + Ус (дар) 
M И К 1 м 
< oe IEn? At + 12117 + 1 эм") 


+C{ 


|ы] олы» 


м м 
(At)? + У "2 + У [т?л 
п=1 n=l 


M M 
+ Уг" 2де У Io" — пФ" |де Н—%(Н* + At)'}. (60) 


n=1 жЕ 
在 上 式 中 ， 令 M = 入， 应 用 离散 的 Gronwall 不 等 式 (15)-(17) 可 以 推 得 
|| < C(At + ak + н?®*-1), 
代入 上 式 ， 再 次 应 用 Gronwall 不 等 式 得 到 


м 
e+ (У даруу) < c(t + ak + HY), (61) 
n=1 


由 关系 式 
lu” — i < n+ 1, E” — ÈRI < lle + E” — Г 

№ (15), (16) 即 可 得 到 定理 的 结论 

定理 3.1 和 定理 3.2 说 明了 两 重 网 格 算法 一 个 明显 的 事实 : 为 了 求解 细 网 格 A 上 的 混合 有 
眼 元 解 ， 可 以 采用 两 重 网 格 算法 ， 选 取 粗 细 网 格 的 步 长 满足 关系 式 h* = OUE2k-1)， 两 重 网 格 
算法 的 解 与 原 非 线性 迭代 算法 的 解 保持 同样 的 计算 精度 。 但 是 ， 由 于 两 重 网 格 算法 只 在 粗 网 格 
上 进行 非 线 性 迭代 运算 ， 在 需要 求解 的 细 网 格 上 只 进行 一 次 线性 运算 ， 这 样 就 保证 了 计算 次 数 
的 减少 、 计 算 效 率 的 提高 。 下 面 的 算 例 将 进一步 验证 这 点 。 
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4 数值 算 例 
考虑 非 线 性 对 流 扩散 方程 的 模型 问题 


ёи +b-Vu—V-(aVu)=f(u), (x,t) єх (0,T], 
u(x,t) = 0, (x,t) € дО х (0,T], (62) 


u(x, 0) = u(x), xEQ, 


其 中 Q = [0,1] x (0, 1], x = (2, y)7, f(u) = -u?+G(z, y,t), G(x, y, t) ш (62) KIRAR u(x, t) = 
tzy(1 — 2)(1 — yje” +! 唯一 确定 。 取 a = 0.001, b= (1,1)7, W (62) 是 一 个 对 流 占 优 的 扩散 方 
程 。 

选取 时 间 步 长 At = 1.25 x 10-5， 计 算 20 层 。 按 照 定 理 3.2 的 结论 ， 当 选取 粗细 网 格 直径 比 
为 h* = O(E2 1) 时 ， 由 定理 3.1 可知 特 征 混合 有 限 元 ( 非 线性 迭代 ) 算法 (т) 的 解 和 特征 混合 
有 限 元 两 重 网 格 算法 (38)、(8) 的 解 具 有 相同 的 计算 精度 。 为 了 保证 起 到 加 速 收 敛 的 目的 ， 有 
限 元 空间 多 项 式 的 次 数 上 至 少 需要 取 2， 在 这 里 我 们 就 取 大 = 2 进行 计算 ， 这 时 的 粗细 网 格 比 
为 h = О(НЗ/?). 

采用 算法 (7) 计算 方 程 的 特征 混合 有 限 元 解 ， 需 要 对 区 域 Q 进行 细 网 格 放 分， 在 这 里 进 
行 全 矩形 单元 划分 ， 如 果 取 细 网 格 步 长 h = 1/32, BEAM Ar EWR 32x 32 = 1024 个 
点 的 近似 值 如 ， 求 解 时 要 用 到 非 线性 迭代 。 

如 果 采 用 本 文 的 两 草 网 格 算法 ， 则 按照 关系 式 及 = О(НЗ/2), FR OBLEKA Н = h3/? я 
1/10 和 矩形 单元 剖 分 ， 得 到 粗 网 格 Ay， 用 非 线 性 迭代 格式 (38) Т ЖШН шн. АЕ 
细 网 格 An 求解 线性 问题 (8) GANERRA n HAE ur л 软 。 表 1 和 表 2 分 别 给 出 了 非 线性 
友 代 解 和 两 重 网 格 解 算法 的 计算 时 间 和 计算 误差 。 


Al: BEACH GO, HL? RE Ж2: AERA un th L RE 
H h uå CPU time h He all CPU time 
5 + 44.0 x 1073 10.6” Б 4.40 x 107° 12.1” 
i эз 5.20 x 107° 165.9” 35 2.50 x 107% 230.0” 


由 表 1、2 可 以 看 出 ， 两 重 网 格 算法 在 提高 计算 效率 方面 非常 显著 ， 计 算 时 间 只 占 原 算法 计 
算 时 间 的 三 分 之 二 ， 而 且 两 重 网 格 解 与 原 算法 解 保持 同样 计算 精度 ， 因 此 用 两 重 网 格 解 作为 问 
题 (62) 的 近似 解 ， 在 不 影响 精度 的 情况 下 ， 有 效 的 提高 了 计算 效率 。 

两 重 网 格 算法 为 求解 非 线性 方程 (包括 非 线性 对 流 扩散 方程 ) 提供 了 提高 计算 效率 的 途径 。 
其 关键 所 在 就 是 只 在 粗 网 格 上 进行 非 线性 迭代 ， 这 样 的 计算 结果 并 不 影响 最 终 解 的 精度 。 两 重 
网 格 算法 与 特征 有 限 元 法 的 结合 不 仅 消 除了 数值 振荡 现象 ， 还 提高 了 非 线性 对 流 扩散 问题 的 计 
算 效 率 ， 是 一 种 值得 推广 的 稳定 、 高 效 算法 。 
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Two-grid Method for the Characteristic Mixed Finite Element 
Approximations of 2D Nonlinear Convection Diffusion Problems 
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Resources and Hydro-electric Engineering, Xi’an University of Technology, Xi’an 710048) 


Abstract: A two-grid method for numerically solving the nonlinear convection-dominated diffusion 
problems is presented. This method includes two steps. First, a nonlinear problem is solved on a coarse 
grid, and then a linear problem generated by the coarse grid solutions is solved on a fine grid. The 
theoretical analysis and numerical example show that the new method is numerical stable and efficient. 
Keywords: convection-diffusion equations; characteristic mixed finite element; two-grid method; con- 
vergence 


